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Abstract 

We study the possibility of spatially non-uniform ground state in (1+1)- 
dimensional models with quartic fermi interactions at finite fermion densities 
by introducing chemical potential fi. We examine the chiral Gross-Neveu 
model and the Cooper pair model as toy models of the chiral symmetry 
breaking and the difermion pair condensates which are presumed to exist 
in QCD. We confirm in the chiral Gross-Neveu model that the ground state 
has a crystalline structure in which the chiral condensate oscillates in space 
with wave number 2fi. Whereas in the Cooper pair model we find that the 
vacuum structure is spatially uniform. Some discussions are given to explain 
this difference. 
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I. INTRODUCTION 



Illuminating the vacuum structure of QCD at high baryonic densities, which should be 
useful to understand physics of compact stars and heavy ion collisions, remains challenge 
for particle and nuclear theorists. One of the most interesting in the related subjects is 
the color superconductivity of quarks |],|2||1- Quark Cooper pairs condense due to the BCS 
instability of Fermi surface at high densities and low temperatures, while chiral condensates 
is dynamically generated by particle-antiparticle pairing at low temperatures and densities. 
Since quarks have a color component of the SU (3) gauge group, the ground state of QCD 
is a superconductor with spontaneous breaking of color symmetry as a result of diquark 
condensates. 

In this paper, we address the question of spatial variations of fermionic condensates, the 
problem first investigated by Deryagin, Grigoriev and Rubakov (DGR) in QCD with infinite 
number of colors N c ||. Working in the perturbative regime g 2 N c <C 1 (g is the gauge 
coupling) and by using the variational method, DGR discovered a previously unsuspected (in 
the context of QCD) instability of the Fermi surface. Namely, they found that an instability 
exists in the large N c QCD at finite densities which triggers condensate of particle-hole pairs 
(Overhauser effect) having almost the same momentum p with modulus \p\ = /i, where 
fi is the chemical potential. They showed that the ground state of the theory is spatially 
nonuniform, and has a periodically varying chiral condensate with wave number 2/x. 

Following DGR, analyses which utilize either renormalization group analysis || or Wilso- 
nian effective action formalism || as well as the Nambu-Gorkov formalism were attempted 
to examine which effect, the Overhauser or the BCS effects wins at large, at finite N c . It was 
shown in these analyses that the BCS effect overtakes the Overhauser effect in weak cou- 
pling regime, while Overhauser wins against BCS in strong coupling regime. More recently, 
the possibility of crystalline ground state was investigated in the context of color supercon- 
ductivity H||[10|. The underlying mechanism which leads to the crystalline structure is a 
mismatch in the fermi surface of up and down quarks, and is apparently quite different from 
one for the periodic chiral condensate discussed by DGR. 

To have better understanding of the relationship between these related but different 
phenomena, it is desirable to have a theoretical laboratory to investigate the underlying 
mechanism of these different phenomena. In this paper, we examine two simple (1+1)- 
dimensional models with quartic fermi interactions which would serve for this purpose. We 
examine the chiral Gross-Neveu model [|ll|] and the Cooper pair model proposed by Chodos, 
Minakata and Cooper |12 . 



The chiral Gross-Neveu model [JTTJ] is a natural choice as the two-dimensional model 
which displays dynamical breaking of continuous chiral symmetry. It is also well known 
that the Gross-Neveu model at infinite number of flavor has a number of similar properties 
with QCD; the renormalizability, the asymptotic freedom, and the similar phase structure 
on fi-T plane with those suspected for the two-flavor QCD. 



The Cooper pair model |0| is a toy model for "color superconductivity" in QCD but 



without color degrees of freedom. It nevertheless serves as a renormalizable asymptotically 
free model field theory which admits Lorentz scalar fermion pair condensate. The model 
possesses O(N) flavor symmetry and U(l) symmetry of fermion number, and the authors of 
Ref. |T!| have shown that at the large- N limit the U(l) symmetry is dynamically broken by 
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the Lorentz singlet di-fermion condensate. It was also shown that the model, when coupled 
with the Gross-Neveu model, has a very similar phase structure with that of QCD with two 
flavor of quarks [[13[ . 

Strictly speaking, continuous symmetry such as the U(l) symmetries of chirality or 
fermion number cannot be violated by the Coleman- Mermin- Wagner theorem JTJJ in (1+1) 
dimensions. However, it is well understood by now in what sense the theorem is "invalidated" 
in the large- N limit. As Witten argued [|15|| , the correlation functions of order parameters 
have power law behavior at large distances, and the power is inversely proportional to 
N. Therefore, if we take the limit iV — > oo first and then send spatial distance between 
two order parameters infinity to use the cluster property, we formally obtain nonvanishing 
vacuum expectation value of the order parameter. The procedure seems to give a consistent 
field theory, which is usuable as a theoretical laboratory for investigating various properties 
of condensate phenomena. 

We develop a systematic method for investigating the possibilities of spatial varing con- 
densates, and then apply the formalism to discuss the chiral Gross-Neveu model and the 



Cooper pair model in a unified way. In a recent paper, Schon and Thies [lb] treated the 
chiral Gross-Neveu model and concluded that the model has chiral wave condensates with a 
wave number 2/z. Our formalism differs from theirs on several respects, in particular on how 
to regulate ultraviolet divergence. We will reexamine the chiral Gross-Neveu model in the 
light of our formalism and confirm their results. On the other hand, We obtain a spatially 
uniform fermion pair condensate as the stable ground state of the Cooper pair model, in 



agreement with JT2 



In Sec. [Tl] we introduce the quasi-particle picture and the Hamiltonian is diagonalized 
in terms of the quasi-particles. In Sec. [I I J the effective potential is derived by taking the 
expectation values of the Hamiltonian with respect to properly defined variational vacua 
and we search for a vacuum solution. In Sec. [TV] we relate the quasi-particles to the free 
particles through the Bogoliubov-Valatin transformation, which leads us to the fact that 
the spatial variations of the chiral condensate yield due to the nonvanishing momentum of 
particle- hole/antiparticle pairs. We also discuss validity of the ultraviolet regularization by 
considering the fermion number of vacuum. In Sec. [V] the similar discussion to the chiral 
Gross-Neveu model is applied to the Cooper pair model. Sec. [\T| is devoted to conclusions. 
In Appendix we discuss the ultraviolet regularization scheme in detail. 



II. THE HAMILTONIAN OF CHIRAL GROSS-NEVEU MODEL 

First of all, we will focus our attention on the chiral Gross-Neveu model in which the 
Lagrangian is given by 



C 



/ + 7 VJ A + ^ 



^(l + 7 5 )^ 



(1) 



Here, ipi is a two-component spinor with a flavor index of SU(N). Repeated flavor indices 
in Eq. ([I]) are meant to be summed. To solve the model, N is sent to an infinite, keeping 
Ao = g 2 N fixed. We also introduce a chemical potential /i which is a constant source for 
fermion number in order to take account of the finite fermion density of the system. The 
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Lagrangian has a U(l) chiral symmetry, which is dynamically broken at zero temperature 
and /i = [[ll|]. Under the assumption of spacially uniform ground state the symmetry is 
shown to be restored at finite temperature and/or chemical potential ||17|| . A recent work 



that relaxes the assumption of spatial uniformity to accomodate a periodically varing one 
fl6|j has shown that a chiral crystalline state with (^(1 + 7 5 )V'i) oc e 2lfix has lower energy 
than uniform chiral symmetry breaking or unbroken state. Using this model, we will first 
develop a method for investigating possibilities of spatially varing condensates, which can 
also be applied to the Cooper pair model. 

Our strategy is to use the auxiliary field method. Following standard techniques [TT 



we 



can add the term involving complex, auxiliary field cr(x) without affecting on the dynamics 
of the theory: 



C -> C- At 

2g 2 



a + ^(l + 7 5 )^ . (2) 



An effective action of the model is obtained by doing the functional integral with respect to 
both auxiliary and fermion fields. In the large- N limit, however, a(x) can be regarded as a 
classical, background field for fermions because it is determined by the method of stationary 
phase. Instead of integrating out the fermion fields in the effective action, we diagonalize 
the fermionic Hamiltonian in the presence of the background field a(x) in terms of the 
creation and annihilation operators of quasi-particles. The free energy is then obtained as 
expectation values of the Hamiltonian with respect to variational ground states for each a. 
At this moment, since a U(l) chiral transformation including the auxiliary field is indicated 
by ipi — > e lOL1 ipi and a — > e 2ia a (a is an arbitrary constant), non- vanishing classical field a 
implies spontaneous symmetry breaking of the U(l) chiral symmetry. 

To analyze the chiral Gross-Neveu model and the Cooper pair model in parallel, it is 
convenient to introduce the charge conjugation field defined by 

r t (x) = c^ T (x), 

where C is a 2 x 2 matrix with the properties (see, for example, Bjorken and Drell fT8|| ) 

C^rfC = -Y T , C = -C T . 

In our convention of Dirac matrices, 7 = G\, 7 1 = — io<i and 7 s = cr 3 (er's are Pauli matrices), 
C = — 7 1 and charge conjugation field is described as 

tf(x)=>ftf(x). (3) 

By using the integral identity 

X - J d 2 x^ + 7 °/i) i>i = \\ - 7 V) C (4) 

the fermionic Hamiltonian density corresponding to the Lagrangian of the chiral Gross-Neveu 
model can be written as 
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dC , dC ,„ 



dipi 



\cr\ 



2<r 



+ 



1 






"ft " 






2 






h c 







(5) 



where the first quantized Hamiltonian H 



h 

h c 



is given by 



h = —ijd x — [i + 7°Re a — ryTm a 
h c = -ry 5 <9 x . + n + 7°Re a + i^lm a 



(6) 



Note that the sign of Ima term flips because the charge conjugation is a kind of chiral 
conjugation (see, Eq. (^)) and the imaginary part of an auxiliary field a transforms as 

Our next task is to identify the eigenstates of Eq. (|5|) and to expand the fermion field 
in terms of them. To do this, we assume that a(x) has the following static standing wave 
form: 



a(x) 



en e 



2iKx 



(7) 



where |<r| and K are arbitary constants which will be determined later. We solve eigenvalue 
problem = u n ^ by Fourier expanding with the chiral phase factor e ±tKxi . We obtain 
a complete, orthonormal set of eigen functions \1/ = [ip n (x), ^(x)] T . It consists of four types 
of continuum state labeled as n = (p, a), (p, b), (—p, a), (— p, b) with the following eigenvalus 
(p is the Fourier mode and we call this type of eigenstates "branch", e.g., a-branch): 



Upa - 
Upb = 
^-pa 
U-pb 



p 1^* 
p ~\~ /^* 



1ppa{%) 

i) pb (x) = 



3 «(p+7 5 K)x 



-pa 



pb 



f) 

x) 



«(p+7 5 K)x 



U_ 



^pb(x) 



7 5 K)x 



pa 



Up, 

e i(p-f 5 K)x u c 



(8) 



pb ' 



0. 



Here, e p = yp 2 + \o~\ and //* = fx — K. The two component spinors in Eqs. (|) are 



u, 



u c _ v = 7 V . 



cos \6p 
sin \Q V 

sin \Q V 

— cos \Q p 



1 


V e p + P 




_ sJtp-V _ 



l 


V e P-p 


V 7 ^ 


_~vtp + p_ 



(9) 



Note that the wave functions are no longer momentum eigenstates in the presense of K 
because of the chiral phase factor. While cj's may be regarded as the energies of a fermion 
or an antifermion with momentum p, mass |cr| and chemical potential /i*, they must be 
interpreted merely as quasi-particle dispersion relations (especially, p is not a momentum of 
particles). 

The fermion field can be expanded in terms of eigenstates (||). The coefficients of four 
branches of eigenstates are not independent with each other, but a (&)- and a(6)-branches 



5 



must be hermitian conjugate with each other to assure that if)f = ^ip*. The second quantized 
fermion field is then given by 





f dp 


e i(p+^K)x( 




~ J 2n 


e i(p-j 5 K)x/ 



U_pCL_pl ~\~ Upbpi j 



(10) 



Here, a P i, b p i and their hermitian conjugates are creation and annihilation operators of quasi- 
particles. To hold the canonical anticommutation relations, {ifji(x),ifjj(y)} = Sij5(x — y), 
they must obey anticommutation relations; {a pi ,a q j} = {b P i,bgj} = 2n5(p — q)5ij, and all 
other anticommutators vanish. 

Using the decomposition into Fourier space, the field operators can be expressed in terms 
of quasi-particle operators, a's and 6's. For example, Hamiltonian and fermion number 
operator become 



2A 
dx 



H = J dxH 

1 f dp / t : 

+ 2 J 2~K \ UJ P a °'P ia P i ^ ^pWpibpi 

~~ 2 J 2ti \ P iapi pi i 



nV- 
i 



+ <^-pa a -pi a -pi + UJ - 



pb b- P ibl pi ) , (11) 



a — pi b—pib—pij ) 



(12) 



p% ^_ pi 

where L — J dx is the spatial volume. It is found that a pi is the annihilation operator of 
the a-branch simultaneously with the creation operator of the a-branch. It should be also 
noticed that the quasi-particle operators diagonalize the fermion number operator differently 
from the Cooper pair model. This becomes the crucial point to distinguish between the two 
models, as one can find below. 



III. EFFECTIVE POTENTIAL AND VACUUM STRUCTURE FOR CHIRAL 

GROSS-NEVEU MODEL 

In this section we want to construct and analyze the effective potential for the chiral 
Gross-Neveu model. We begin by calculating expectation values of the Hamiltonian with 
respect to variational ground state which is defined by filling the negative energy states; 

u pa = -^pa > : a pi \a) = , u pb = -u pb > : b pi \a) = 
u pa = -uj pa < : ajjja) = , u pb = -u pb < : b ] pi \a) = ' 

Evaluating the expectation value of the Hamiltonian with respect to \a), the effective po- 
tential becomes (for brevity, we denote (O) = (a\0\cr)) 



V(\a\,K)= {H) 



NL 
| a | 

2^ 2 ./ 2;-;" '■" ' 2 J 2it~- po ' J 2~~ / "' ' J 2tt 

U)pa<0 ^pb<0 



1 r dp 1 r dp r dp r dp 

+ n / T^U-pa+n / 7T-U-pb+ / ^pa + / ^ pb , (14) 
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where the region of integration in the second and the third terms in the last line consists 
of whole range of p. The last two integrals come from "no-pairing" regions in which the 
modes does not give any contribution to dynamical generation of chiral condensate (see, e.g, 
Eq. (0)). Under the condition for absence of no-paring regions, /i 2 < |cr| 2 , there is no effect 
from free fermions and these integrals vanish. 

Eq. (|14D contains several divergences at ultraviolet regions; for example, if we introduce 
the naive Fourier mode cutoff (as we took care before, p is not a momentum of quasi-particle), 
\p\ < A/2, then the divergent terms come from first two integrals in Eq. ([14]): 

A/2 2 2 

dp A u* , en , A \a\ ... 

fu_ pa = -— + - U- In — y - JJ- + O 1/A 

Z7T 671 Z7T 47T (J 27X 

dp A 2 , |cr| 2 , A 2 I cr 1 2 ... , , 

■£ u -* = -^-£ A -^r ln ^-^r + ° (1/A) (15) 

Z7T S7T Z7T 47T (J 

-A/2 1 1 

The quadratic and linear divergences depend on the cutoff scheme, although the logarithmic 
divergence does not. These divergences therefore lead to nontrivial issue for the ultraviolet 
regularization. We will show at the next section that the naive Fourier mode cutoff has 
problems with the number of fermions floating in the vacuum. We employ a different cutoff 
scheme which is free from the trouble. As will be explained somewhat in detail in Appendix, 
it involves the branch dependent cutoff prescription: 



-A/2 
A/2 



afaVbranch : 



2 '2 



6(6)-branch : 



±-k£+k 

2 ' 2 . 



(16) 



The reader may wonder why one should take such branch-dependent shift on the ultraviolet 



cutoff. In fact, it is shown in Refs. |19[ that the similar shift resulting from a phase shift in 
the normal modes of fluctuations around the soliton gives correct quantum corrections to 
the kink mass in the two-dimentional 4 -theory. 

By using the cutoff prescription (|16"D , first two integrals in the effective potential flMQ 
become 

2ti 4tt (j 2tt 2tt 2tx 

-A/2+K 



dp 


A 2 


Z7T 


8tt 


dp 


A 2 


—Co >_ p l = 
Z7T 


8tt 



A / 2 + K , a9 i ,2 a 9 I |2 9 9 

' ''" U ^ A H , A 2 |cr| 2 , /i 2 /I 2 



-A/2-K 



2.- 4. ln ^-17 + g-& + ( 1 / A ) < 17 » 



One notices by comparing Eqs. ( |TTD to ( [3~5|) that additional finite contributions arise from the 
lineary or quadraticaly divergent terms due to the shifts in cutoff. The quadratic divergence 
in Eq. ( O ) is an uninteresting constant which can be subtracted off and the linear diver- 
gences are exactly cancelled between two integrals. The remaining logarithmic divergence in 
Eqs. ([T7|) can be removed, as usual, by the renormalization of the coupling constant. Since 
they contain only |er| dependence, it is enough to carry it out in the case of K = n = 0. 
Therefore, we renormalize by demanding that the renormalized coupling constant A satisfies 
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d 2 V 

9kF 



| a | — m,Q 

K=ft=0 



1 
A 



(18) 



Here \a\ = m designates an arbitary renormalization point on which the coupling constant 
depends. Using this conditions to solve for Ao in terms of A yields the renormalized form of 
the effective potential: 



2 ^ 1 |2 


: V = 


1 |2 
(J 

47T 


"2tt 

_T 


1 I 2 " 

- 3 + In 

mg 


2vr 


2tt' 


"1 


2 ^ 1 |2 


: V = 


47T 


[ 2tt 

l T 


-3 + ln4 


+ \ 


1 |2 

1 



(19) 



+ 



K 

2tt 



1 



2tt' 



When = 0, Eqs. (|19"D turns into the well-known effective potential assuming the homo- 
geneity of chiral condensates |T7|]. If we set K — fi, the effective potential is reduced to one 



in Ref. 16 



It is easy to show that \a\ 



m e 



l-w/X 



is the solution of the gap 



equations derived from the effective potential, Jj^y = ^ = 0. It give the minmum value of 
effective potential as one can see in Fig. [I]. Therefore, the vacuum structure is characterized 
by 



a {x) = (m e l -* /x ) 



(20) 



Using the field decomposition (|T0l) , one can evaluate the chiral condensate of vacuum ex- 
plicitly, 



(^(1+7 5 M> 



-N \a\ e 2iKx — 
1 1 2tt 



dp 



dp 



U pa <0 



W p b<0 



dp 

€ 



V J 



= -Na{x)/\ , (21) 

where the vacuum value of a(x), Eq. (^), is used in the second line of the equation. The 
equation guarantees to hold the equation of motion for cr(x) derived from the Lagrangian 
(|). Consequently, we have reproduced the result in Ref. Jig which the phase standing chiral 
density wave with a wave number 2{i is stable as compared to uniformally chiral symmetry 
broken phase or the symmetric phase. Moreover, we have demonstrated that the condensate 
of the wave number 2/j, arises as the unique minimum of the effective potential. We present 
arguments to clarify the reasons why the wave number is 2/j, at the next section. 



IV. RELATION TO THE OVERHAUSER EFFECT 

To understand why the chiral condensate oscillates in space with wave number 2/x, we 
start by considering that the field expansion, Eq. fllHD, has a connection with an interaction- 
free description through the Bogoliubov-Valatin transformation. For brevity, we drop the 
flavor indices and set L to unity. Setting \a\ = K = in Eq. fllHD, the free field expansion is 
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V>0) = / — e 



a(°)e(g) + 6 (0 ] t e(-g) 
a (o)e(-g) - feife(g) 



(22) 



where the subscript (0) implies the operators corresponding to free massless fermion and 
antifermion, and 6 denotes the Heaviside step function. By equating Eq. fl22|) to Eq. flTOD, 
the Bogoliubov-Valatin transformation reads 

- ™ c 1 f) A(°)t • In 
*~p - LI p+K COb 2 C p U +K bill 2° pi 

bp = b p °l K cos \9 P - a^_ K sin \0 P , 

(0) i/j . (0) -l/i 

- Up+K C0S 2 P a p-Jf Sm 2 P' fOQt 

->)t ros I/9 _ a (°)t sin Ifl lZdj 

- cob 2 u p u_p_ K bin 2 c p , 

a p = b^p_ K cos |6>p + Oyi K sin |0 p , 
fe P = a^J +K cos ±0 P + bfl K sin |0 p . 




p < -K : 



We can then relate the variational ground state |er) to the normal state |0). The latter is 
annihilated by (a^ 1 ") and b^ above (below) the Fermi surface: 

\q\>p: af\0) = bf\0) = 0, 

\q\<p : aft|0) =6(°)|0) = 0, 1 } 

(We implicitly assume that p > such that the fermionic Fermi surface is formed). For the 
case of K = p in particular, one is able to show that |er) can be written in a form of the 
BCS state, 

i*>= n v n n ^ +t n^ +t io), (25) 

P<-H -n<p<0 0<p<^t M<p 



where 



— cos 2 6 p + sin 2 ^p a p+fJ ) -p\p,j 



,(°) „(°)t 

(26) 

-»p 01J.1 2 i/p t ^wo 2 K/ p Uj p+n u 'p-fi- 

Ap 1 " = sin ±0 p + cos ~0 p &r^4-!«> 



Bp = cos ^p + sin ^Opajpl^ap 
= sin i0_ + cos ij0 o aj2./i£S.. 



It can be easily verified that |cr) in Eq. ( Pq) satisfies the definition of it, Eq. (|T3|), by using 
Eqs. ( P^D and (P^D- In the region of p > //, |er) comprises with the particle-antiparticle 
pairing as in the usual chiral condensate vacua but with pair momentum 2p in our case. On 
the other hand, in the region of < p < p, \a) comprises the particle- hole pairing with the 
momentum (jp + p) and — (p — p), respectively. The manner of pairing is depicted in Fig. |2|. 
This figure tells us that either the particle-hole or the particle-antiparticle pairing across 
the Fermi surface replaces particle-antiparticle pairing across the Dirac sea that generates 
homogeneous chiral condensates at p = 0. Exciting the pair of a particle with momentum 
p and a hole with momentum p (p = 0) costs zero energy. Since the pairing wave function 
tan 1 6p at p > has the peak at p = 0, most likely, the zero energy pairing is generated 
when the interaction is truned on. This is the Overhauser instability first proposed for spin 
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density wave of nonrelativistic electron system ||20|| , and recently revived for chiral density 
wave in QCD ||,[5]|||7]]. Because of effects of this instability, the vacuum must have the chiral 
crystalline structure with a period 7r//i. 

Eq. ([25]) tells us an another important result. Since |0) has nonzero fermion number 
owing to the existence of the Fermi sea and the pairing preserves fermion number, |a) must 
also have nonzero fermion number. The expectation value of a fermion number operator 
becomes 



A/2-K A/2+K 
-A/2+K -A/2-K t*V<0 



2vr 



uj pb <0 



2tt 



nl'1 

IT 

7T 7T \ 



| (7 1 



for /x* < \ a \ 



for /i^ > | cr | 



(27) 



by our cutoff scheme (|lq) . At /x* = 0, or equivalently K = /x, it yields the correct fermion 
number of the state \cr), (J°) = NL/jl/tt. We emphasize that the consistent result comes out 
because of the cutoff procedure ([lj]) we employ. Vanishing fermion number would result if 
the naive Fourier mode cutoff were used. 

We comment out about the relation of Eq. ( p7D to early works relating on the zero 
chemical potential system, finally. In the chiral sigma model in which a(x) turns into the 



true backgroun field, the authors of Ref. [21 has shown that the vacuum expectation value 
of fermion number density is = — , where a(x) = \a\ e~ tip( - x ^' y5 , by using the adiabatic 
method. To compare with our model we set tp(x) = —2Kx, so that the fermion number 
density becomes the first case of Eq. (p7|) expecting the factor iV coming from the number 
of flavor. Since the adiabatic method implies the slowly varing condition (d x ip -C \cr\), the 
secound case of Eq. (|2~?D is excluded from this condition. 



V. THE COOPER PAIR MODEL 



We discuss the Cooper pair model in this section, replacing the interaction term in Eq. (IT]) 

2 



with \g 2 



Mi 



c 



A (i 



+ 7V 



+ 2 9 



(28) 



UJI, due to our convention of Dirac 



The Lagrangian fl28|) is the same as the one in Ref. 
matrices. The Lagrangian has a U(l) symmetry which implies the conserved fermion number 
in the model. As was shown in Ref. |12j under the assumption of spatilly homogenuity, a 
Lorentz scalar Cooper pair condensate is dynamically generated in this model; (ipiipi) ^ 
(ipf transform as ipi under the Lorentz transformations). 
We introduce an auxiliary field A(x) by 



2<T 



A + g 2 ^ t 



(29) 
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When one thinks about the Cooper pair condensate, it is convenient to regard the charge 
conjugation of ipi as independent variable of the theory (see, for example, Refs. |l|j2^|). By 
introducing the kinetic term of charge conjugation field by Eq. (|]), we obtain the Hamilto- 
nian density in the following form, 



n 



|A|' 
2? 



+ 



-i>y 5 d x - fi 7°A 
7° A* -vfd* + fj, 



(30) 



As in the case of chiral Gross-Neveu model, we impose the standing wave ansatz for the 
background (auxiliary) field: 



A(x) 



IAI 



2iKx 



(31) 



A(x) is equal to the Cooper pair condensate {tpfipi) apart from a constant factor, so that 
nonvanishing |A| indicates the symmetry breaking of fermion number £7(1) symmetry. 

Eigenvalues of the first quantized Hamiltonian in Eq. ( |30"D consist of the following four 
branches: 



-P+ 



K 



-pa 
pb 



-6 p _ - K, 

-e p+ + K, 



(32) 



where e p = yp 2 + |A| is a energy of a particle with momentum p and mass |A|, and p± = 
p =F fi are momenta relative to the Fermi surface. We can then diagonalize the Hamiltonian, 
and fermionic field is expanded in terms of the eigenstates with phase factor, e ±lKx : 



dp 
— < 

2vr 



,i(p+K)x 



dpi cos \6 p+ + bLpi sin |6> p+ 

Qj COS ~^0p_ 



d P Jh 



i(p—K)x 



2tt 



bpi sin \9 P _ 

bpi cos \0 V _ + a)_ pi sin \Q V _ 
dpi sin 2^p+ 



b-pi cos \0 P+ 



(33) 



The wave functions cos \Q V and sin \Q V are defined analogously in Eq. (|9|). The Hamiltonian 
is obtained in exactly the same form as in Eq. 



TT|) with the energies of Eq. 



Decomposing the fermion number and current operator in the basis of Eq. (p3|) manifests 
the difference from the chiral Gross-Neveu model. In the Cooper pair model, fermion current 
is diagonalized by the quasi-particle operators as the following equation, but fermion number 
is not; 



J 1 



dx 



2 J 2tt ^ 



Qjpidpi 



u P t u pi 



-pv^-pi 



-pi u -pi 



(34) 



The cutoff scheme is the same one as the chiral Gross-Neveu model, because the field 
expansion ( p3|) must reduce to the free theory even if | A| = ji = and K ^ (see Appendix). 
Variational ground state |A) is defined by filling the negative energy states in the same way 
as Eq. (|i~3"D, and then the effective potential becomes 
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V(|A|,AT;//) 



<#> 
NL 

2 k/2-K k/2+K 

|A| 1 r dp 1 r dp 

2\ Q + 2 J 2^ + 2 7 2^ 

-A/2+K -A/2-K 



+ 



dp 



u p b<0 



dp 
2^ 



(35) 



In (^), (if) implies to take expectation value with respect to the state |A). Performing the 
integrations explicitly and using the renormalization condition ( |i~8D (replacing |<r| to |A] in 
the equation), we obtain the renormalized effective potential as 



K 2 < |Ar : V 



K 2 > IAT : V 



A| 2 


"2tt 


-3 + ln 


A| 2 








47T 


_T 


ml 


A|^ 


I 2tt 


-3 + ln 


K 2 








All 


1 T 


ml 



2tt 



2tt' 




1 - 



K 2 



(36) 



K 2 
+ 2^ 



|A£ 
iT 2 



if! 

2tt' 



It is easy to observe that the effective potential is identical with that of the chiral Gross- 
Neveu model, Eq. (0), under the replacement of /i* by K. Then, the vacuum is obtained 
at K = 0, and the Cooper pair condensate of vacuum is given by 



A{x) 



m e 



1-tt/A 



(37) 



(The last equality is straightforwardly proved, similarly to Eq. (^T|) in the chiral Gross-Neveu 
model). Therefore, the ground state with spatially uniform Cooper pair condensate obtained 
in Ref . 



gives the true vacuum in this model. 
The question we have to ask here is why the Cooper pair condensate is uniform differently 
from the chiral condensate in the chiral Gross-Neveu model. In order to answer this question, 
we had to better examine the fermion number and current of the variational ground state 
| A). After some calculations, they turn out to be 



(J°) = NL 



7T 



NL- 



K 



, K K |Ar 



for K 2 < |A| 2 
for K 2 > |A| 2 



(38) 
(39) 



The fermion number is independent of the parameter indicating the spatial varidity K in 
constrast to the chiral Gross-Neveu model. The spatial homogeneity can be explicitly broken 
in the chiral Gross-Neveu model because the fermion number is an order parameter for this 
symmetry breaking and it has a source, namely chemical potential. In the Cooper pair model 
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this symmetry breakdown, if it occurs, must be caused spontaneously because the fermion 
number is no longer the order parameter. Then, the spontaneous symmetry breaking of the 
spatial homogeneity does not occur in our models. 



VI. CONCLUSIONS 

In this paper we have analyzed the vacuum structure of (l+l)-dimensional models with 
quartic fermi interactions at finite densities, with particular emphasis on the possibility of 
spatially nonuniform ground state. Our treatment is restricted into the case of large- N limit 
and of zero temperature. This is the unique place where we can have nonvanishing order 
parameter in 1+1 dimensions in the sense explained in Introduction. 

We first examined the chiral Gross-Neveu model. We have constructed the effective 
potential by allowing the chiral condensate to vary periodically in space. The quasi-particle 
basis diagonalizing the fermionic Hamiltonian of the theory was obtained by the Fourier 
expansion after performing the chiral rotation. We have constructed an effective potential 
and have shown that the vacuum has a crystalline structure in which the chiral condensate 
oscillates with the wave number 2/x. This is consistent with the result obtained by Schon 
and Thies (ITS . 



In the course of constructing the effective potential, we have found that the branch- 
dependent ultraviolet cutoff procedure is required to derive the correct effective potential. 
In the cutoff scheme, the shifts from the explicit cutoff parameter play important role for 
the purpose of providing the system with the correct fermion number. While the branch- 
dependent ultraviolet cutoff procedure may look ugly, it is likely that it is required to preserve 
the fermion number in the system by cancelling the spectral flow induced by chiral rotation 
performed to obtain the basis used in our treatment. 

We have also constructed the variational ground state in terms of free particle basis and 
performed the Bogoliubov-Valatin transformation to obtain quasi-particle basis which diag- 
onalize the Hamiltonian. We have shown that the wave number 2/i of the chiral condensate 
results from the Overhauser instability which generates particle-hole/antiparticle pairs with 
the total momentum 2/i, or —2/i. 

We have examined the possibility of spatially varying ground state in the Cooper pair 
model proposed by Chodos, et al. |I2] by allowing the Cooper pair condensate to vary 
periodically in space. By constructing and analyzing the effective potential, we have shown 
that the ground state with spatially uniform Cooper pair condensate obtained in Ref. [O 
gives the true ground state of this model. 

We should note that our model, which contains a single chemical potential common to 
iV-flavor fermions, does not quite mimic the QCD motivated (3+l)-dimensional model in 
Refs. |||9,10| in which the mismatch in the Fermi surfaces of up and down quarkes are 



the essential ingredient for the crystalline color superconducting phase. It would be very 
interesting to construct a (l+l)-dimensional model which possesses the analogous structure 
and work it out in the similar way as done in this paper. However, the formulation of the 
large- N field theory with two different chemical potentials is highly nontrivial and we did 
not enter into the problem in this paper. 

The spatial variations of Cooper pair condensate, if it forms, is due to the condensation of 
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particle-particle pairs (or antiparticle-antiparticle pairs) with nonvanishing total momenta, 
as we observe from the Bogoliubov-Valatin transformation. It means that there is nozero 
expectation value of fermionic current. In other words, the fermionic current becomes an 
order parameter for the breakdown of the spatial homogeneity. 

On the other hand, the role of order parameter is played by fermion number in the chiral 
Gross-Neveu model, so that the breakdown of the spatial homogeneity is caused explicitly 
by the chemical potential (in fact, the symmetry is restored at ft = 0). In the Cooper pair 
model, however, this symmetry breakdown is caused spontaneously if it occurs. The result, 
therefore, shows that the spatial homogeneity is not sopntaneously broken in our models. 
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APPENDIX A: ULTRAVIOLET REGULARIZATION 

In this Appendix we explain the ultraviolet cutoff procedure given in Eq. (|16|) in detail. 

If we set |er| = in Eq. (|7|) in the chiral Gross-Neveu model, it becomes a free field 
theory. In this case the quasi-particle description should reduce to the free field description. 
Eq. (|TUD, however, differs from Eq. fl2"2|) by the factor e ±tKxi . We will show below that the 
superficial difference between the quasi-particle and the free-particle description is removed 
by the ultraviolet cutoff prescription. 

From the Bogoliubov transformation (|23"D at |cr| = we first obtain 



<h = a p °l K for p > 0, 

(0) f ^ n v A1 J 

a p = a p _ K tor p < 0. 

It is natural that the free field description possesses naive momentum cutoff, \q\ < A/2 
in Eq. ([22"D. We observe from Eq. ( |Al|) that the cutoff in the a- branch becomes —A/2 < 
p — K = q < —K and K<q = p + K< A/2, i.e. \p\ < A/2 — K. The remaining mode of 
free fermion, —K < q < K, comes from the 6-branch: 



b p = a m p l K for < p < K, 
bp = -a { °ll K for -K < p < 



(0)f r , ( A2 ) 



Finally, the rest of the Bogoliubov transformation, 

b p = bf_ K for p>K, 
b p = b p °l K for p < —K, 



(A3) 



leads to the cutoff in the 6-branch, —A/2 < p + K = q and q = p — K < A/2, namely 
\p\ < A/2 + K. In addition, we can verify by using Eqs. (|A"T|)-([A"3"D that the quasi-particle 
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decomposition fllPf) with the cutoff (|16D is reduced to the free-particle decomposition 
with the naive momentum cutoff, |g| < A/2. 

In the above discussion, we have proved that the ultraviolet cutoff must be Eq. ([16]) for 
| a | = at least. We believe that this cutoff is correct one at |er| 7^ because the cutoff 
procedure must be insensitive to the vacuum expectation value, as mentioned in Sec. [TTH ; 
thus, this cutoff is required to provide the state \a) with the correct fermion number. 

Now, we comment on the cutoff prescription from a different point of view. We enclose 
the system in a finite box of length L for the purpose of counting the number of modes 
correctly. We shall impose periodic boundary conditions on ipi(x) and ip^{x). This leads to 
the following conditions on K and p because of the chiral phase factor in wave functions (^j), 

K = ^ P-^U + Ib), (A4) 



L L V 2 

where B and n are integer. One can see from the cutoff fll6|) that n take (2N A — B + l) values 
from —A/a to N A —B for a-branch, and (2N A +B+1) values from —N A —B to N A for o-branch, 
respectively, for fixed B (i.e., fixed K). Here, N A = AL/4ir is an explicit cutoff parameter. 
The number of modes included in divergent sums to evaluate the effective potential ( |i4|) is 
kept constant, independently of B; (2N A -B + l) + (2N A + B + 1) = 2{2N A + 1). We depict 
the cutoff prescription in Fig. |3|. This figure tells us that the cutoff implies summing modes 
in order from above with respect to the absolute value of energy, keeping the total number 
of modes fixed when \x = 0. 

We also note that the fermion number of the state \a) results from the difference of mode 
between a- and o-branch when we take the cutoff as Fig. ^|. Since a-branch carries fermion 
number — | and o-branch carries +| as one can find from Eq. (fP2"D, the fermion number of 
\a) becomes -|(2A/ A - B + 1) + ±(2N A + B + 1) = B times N of flavor, namely, AXf . 
It is consistent with the first equation of Eq. ( |2"TD because we now ignore the region where 

For the Cooper pair model, we can discuss the cutoff procedure as the chiral Gross-Neveu 
model. By comparing Eq. (p3|) to Eq. (p2|), we obtain the Bogoliubov transformation 



( _ (°) lfl _ (°)t • la 

K ^ r. ■ J ap ~ a P+K COS 2°P+ a -p+K Sm 2 P P+' 
P \ u la t(0)f • in 

{ bp = b p l K cos 1 6 p _ - b\' p '_ K sm ^9 P _ , 

-A < £> < A . < (0)f (0)f . (AO J 

I Op = a_pVA' cos \vp- ~ °-p-K sm 2&P- , 

p< _ K . j a p = ° {0) p-k cos \0 P+ + 6p ^ sin |0 p+ , 
' I b p = a_ p + K cos |6 l p_ + a p + K sin |#p_ . 

The condition in which the quasi-particle description is reduced to the interaction-free de- 
scription at |A| = in spite of the presence of K, leads to the cutoff fllfip. 

We also depict the same figure as Fig. |3| because the dispersion relations of quasi-particles 
in the two models coinside with each other at \i = 0. In the Cooper pair model, however, the 
quasi-particle does not carry definite fermion number, but the fermion current. Therefore, 
the fermion current of the state |A) is generated as we vary the value of K, as indicated 
in Eq. (BW). This means from the viewpoint of the Bogoliubov transformation (|A5|) that 
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I A) is composed of two (free)fremions pairing with total momentum IK or two antifermions 
pairing with total momentum —2K. 
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FIGURES 



V 




FIG. 1. The effective effective potemtial for the chiral Gross-Neveu model. Each axis is scaled 
by m = moe 1_7r / A . The contour of the effective potential is depicted on the bottom of 3D plot. 
For the Cooper pair model, replace |<r| and ^* by |A| and K, respectively. 




P 



FIG. 2. (a) The manner of pairing for the chiral Gross-Neveu vacuum. u q is free particle 
energy; thus, ujq J is of the fermion, co^ is of the hole and co^ is of the antifermion. The pairing 
with total momentum 2\i (p > 0) or — 2fi (p < 0) yields, (b) Pairing wave function <p(p). 
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